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ABSTRACT 

To study the fragmentation and gravitational collapse of dense cores in infrared dark clouds (IRDCs), 
we have obtained submillimeter continuum and spectral line data as well as multiple inversion tran¬ 
sitions of NH3 and H2O maser data of four massive clumps in an IRDC G28.53—0.25. Combining 
single dish and interferometer NH3 data, we derive the rotation temperature of G28.53. We identity 
12 dense cores at 0.1 pc scale based on submillimeter continuum, and obtain their physical properties 
using NH 3 and continuum data. By comparing the Jeans masses of cores with the core masses, we 
find that turbulent pressure is important in supporting the gas when 1 pc scale clumps fragment into 
0.1 pc scale cores. All cores have a virial parameter smaller than 1 assuming a inverse squared radial 
density profile, suggesting they are gravitationally bound, and the three most promising star forming 
cores have a virial parameter smaller than 1 even taking magnetic held into account. We also asso¬ 
ciate the cores with star formation activities revealed by outflows, masers, or infrared sources. Unlike 
what previous studies suggested, MM1 turns out to harbor a few star forming cores and is likely a 
progenitor of high-mass star cluster. MM5 is intermediate while MM7/8 are quiescent in terms of star 
formation, but they also harbor gravitationally bound dense cores and have the potential of forming 
stars as in MM1. 

Subject headings: stars: formation — ISM: molecules 


1. INTRODUCTION 

Massive infrared dark clouds (IRDCs) in the Galaxy 
have been recognized as the birthplace of high-mass 
(M > 8 Mq) stars, since they are massive, cold, and 
dense (Simon et al. 2006; Pillai et al. 2006; Rathborne 
et al. 2006). Millimeter/sub-millimeter interferometric 
observations toward IRDCs have found 0.1 pc scale frag¬ 
ments or cores that might harbor high-mass stars (Zhang 
et al. 2009; Zhang & Wang 2011; Wang et al. 2011, 2014; 
Beuther et al. 2013; Peretto et al. 2013). These cores 
are often associated with molecular outflows, but lack 
complex organic molecular line emission that are tracing 
hot cores. Therefore, they are the most promising speci¬ 
men for studying the very early phase of high-mass star 
formation. 

One important question concerning these cores in 
IRDCs is how they fragment and evolve, until high-mass 
stars are born. Zhang et al. (2009) found that in IRDCs 
parsec-scale clumps fragment into 0.1 pc scale cores, and 
Wang et al. (2011, 2014) found that 0.1 pc scale cores 
further fragment into <0.1 pc scale condensations, all 
of which are more massive than the corresponding ther¬ 
mal Jeans masses but are consistent with turbulent Jeans 
masses. Therefore, in IRDCs turbulent pressure is es¬ 
sential in supporting the fragmentation, so that massive 
cores are able to form and grow. Pillai et al. (2011) 
studied 0.1 pc scale cores in two high-mass star forming 
regions, and found that all the cores are gravitationally 
bound, with their gravitational mass far exceeding the 
virial mass, therefore turbulence is not sufficient to stop 
cores from fast collapsing. 

To study the impact of turbulence in the initial frag¬ 


mentation and core growth in IRDCs, we select four 
clumps in an IRDC G28.53—0.25 (G28.53 hereafter). 
G28.53 has a kinematic distance of 5.4 kpc (Rathborne 
et al. 2006) and a luminosity of ~3500 L 0 (Rathborne 
et al. 2010). It has been mapped in 1.2 mm continuum 
with the IRAM 30m single-dish telescope (Rathborne 
et al. 2006, 2010, also see Figure 1), which reveals a total 
mass of ~10 4 M 0 . Ten continuum peaks are identified, 
each of which is a few hundreds of M 0 and of ~1 pc 
scales, therefore they are typical clumps that form mas¬ 
sive stars. Rathborne et al. (2010) also classified these 
clumps as active, intermediate, or quiescent in terms of 
star formation based on the infrared emission. Among 
them, the most massive (~10 3 M 0 ) one at the center 
of G28.53, MM1, is classified as quiescent. A less mas¬ 
sive but more luminous clump at the southern end of the 
cloud, MM5, is classified as intermediate. Two clumps 
in the southeast corner, MM7 and MM8, are classified as 
quiescent. Sanhueza et al. (2012) obtained 3 mm spectral 
lines of the ten clumps and confirmed the classification 
of Rathborne et al. (2010) based on chemistry. 

Interferometric observations at angular resolutions 
of l"-2" toward MM1 revealed further fragmentation 
(Rathborne et al. 2008; Swift 2009). Rathborne et al. 
(2008) resolved a 0.1 pc scale core in MM1 into three 
smaller condensations. Although these condensations do 
not present detectable spectral line emission at an arcsec- 
ond resolution, the single-dish observations of MM1 sug¬ 
gested active star formation given broad line wings and 
detection of SiO emission which is a typical shock tracer. 
Swift (2009) mapped MM1 with the Submillimeter Array 
(SMA) in 345 GHz and found a massive (~60 M 0 ) core 
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known as ‘core 2 ’, which is different with the one found 
by Rathborne et al. (2008). No CO outflows or hot core 
tracers were found in this core. 

However, given its large amount of gas and its loca¬ 
tion at the center of the gravitational potential well of 
G28.53, it is puzzling that MM1 is more quiescent than 
MM5. Previous single-dish observations might miss the 
deeply embedded cores in MM1 that are actively form¬ 
ing stars. The two interferometric studies detected a few 
dust cores but lacked molecular lines to trace the star for¬ 
mation, and the CO emission at 345 GHz as an outflow 
tracer was filtered by the interferometer. Therefore, it 
is worth revisiting the 0.1 pc scale fragmentation within 
MM1, and compare it with the ‘intermediate’ MM5 and 
‘quiescent’ MM7 and MM 8 . 

In this paper, we will use high angular resolution NH 3 
line observations to trace kinematics and determine gas 
temperatures. NH 3 itself is also useful in tracing dense 
gas with a critical density of ~10 4 cm -3 (Swade 1989). 
We will also use millimeter continuum and spectral lines 
to resolve 0.1 pc scale cores in the four clumps and 
trace star formation. The combination of kinematics 
(e.g. linewidth, velocity components along the line of 
sight) and temperatures from NH 3 lines and dust emis¬ 
sion from millimeter continuum will allow us to deter¬ 
mine the properties of cores reliably. We will introduce 
details of the observations in Section 2. Then we present 
the results in Section 3, and discuss the fragmentation 
and star formation in the clumps in Section 4. In Sec¬ 
tion 5 we summarize the main results and present our 
conclusions. 

2. OBSERVATIONS 
2 . 1 . SMA observations 

2.1.1. 230 GHz observations of G28.53 MM5, MM7, & 

MM8 

In 2010, we observed G28.53 MM5, MM7, and MM 8 
with seven antennas of the Submillimeter Array 6 (SMA; 
Ho et al. 2004) in the compact configuration at 230 GHz 
band. Frequency-dependent bandpass solutions were ob¬ 
tained by observing quasars 3C454.3. Time-dependent 
gain solutions were obtained by observing the quasar 
1911—201 every 20 minutes. The flux calibration was 
done using Titan and Neptune. The correlator covers 
rest frequencies of 216.7-220.7 GHz in the lower side 
band, and 228.7-232.7 GHz in the upper sideband, with 
a uniform channel width of 0.812 MHz, which is equiv¬ 
alent to 1.1 kms -1 at 230 GHz. System temperatures 
were 100 120 K and opacity at 225 GHz was 0.06-0.1. 
The observations are summarized in Table 1. 

The visibility data were calibrated using MIR 7 . Cal¬ 
ibrated data were then inspected and imaged using 
MIRIAD 8 (Sault et al. 1995) and CASA 9 (McMullin 
et al. 2007). Continuum emission was extracted by av¬ 
eraging line free channels in the visibility domain, then 
was imaged using combined data from both sidebands. 

6 The Submillimeter Array is a joint project between the Smith¬ 
sonian Astrophysical Observatory and the Academia Sinica Insti¬ 
tute of Astronomy and Astrophysics and is funded by the Smith¬ 
sonian Institution and the Academia Sinina. 

7 https : //ww. cf a.harvard. edu/~ cqi/mircook.html 

8 http://www.cfa.harvard.edu/sma/miriad/ 

9 http://casa.nrao.edu/ 


The spectral lines were cut out from the continuum- 
subtracted visibility data and were imaged separately. 
All images were CLEANed with a robust parameter of 
0.5 to obtain a balance between good sensitivity and side 
lobe suppression. Image properties are listed in Table 2. 

2.1.2. Archival 230 GHz data of G28.53 MM1 

We used the SMA 230 GHz archival data of G28.53 
MM1. The observations were carried out in the 230 GHz 
band in two tracks in 2009. The pointing center is 
at (18 /l 44 m 18.0 s ,—3°59'23.0"). The SMA was in the 

compact configuration with seven antennas in the ar¬ 
ray. Frequency-dependent bandpass solutions were ob¬ 
tained by observing the quasar 3C273. Time-dependent 
gain solutions were obtained by observing the quasar 
1751+096 every 28.5 minutes. The flux calibration was 
done using Callisto. The two tracks both cover 2 GHz 
in each sideband of the SMA, with a uniform channel 
width of 0.812 MHz, which is equivalent to 1.1 kms -1 at 
230 GHz. The tracking frequencies of the two tracks are 
different, so that J = 2-1 lines of CO isotopologues and 
the CH 3 OH line at 229.759 GHz were observed in one 
track and the H 2 CO line at 225.698 GHz was observed 
in the other. System temperatures were ~80-160 K and 
opacity at 225 GHz was ^0.06-0.08. The observations 
are summarized in Table 1 . The visibility data were cali¬ 
brated and imaged in the same manner as in the previous 
section. Image properties are listed in Table 2. 

2 .2. VLA and GBT NH 3 observations 

We observed five positions in G28.53 using the Na¬ 
tional Radio Astronomy Observatory (NRAO ) 10 Very 
Large Array (VLA) in the D configuration in two obser¬ 
vation runs in 2010 January, to obtain the NH 3 (J, K ) = 
(1,1) and (2,2) transitions. The 3.125 MHz bands were 
each split into 128 channels with a channel spacing of 
24.4 kHz, equivalent to 0.3 kms -1 . The two bands si¬ 
multaneously covered the main and two inner satellite 
transitions of the NH 3 (1,1) line, as well as the main and 
one inner satellite transitions of the NH 3 (2,2) line. Co¬ 
ordinates of these positions are listed in Table 1. 

The quasars 0137+331 (3C48) and 1331+305 (3C286) 
were used for flux calibration. Bandpass calibration 
was done with observations of 1229+020 (3C273) and 
2253+161 (3C454.3). Gain calibration was performed 
with periodical observations of 1851+005 every 15 min¬ 
utes. The visibility data were calibrated using CASA. 

The NRAO Green Bank Telescope (GBT) was used to 
observe a 8 ' x 7.5' region in G28.53 in 2010 February, 
to obtain the NH 3 ( J,K ) = (1,1) and (2,2) transitions. 
The data were calibrated and imaged using GBTIDL 11 . 
The NH 3 images were then Fourier transformed into the 
visibility domain based on the VLA visibility models. 

Then the VLA and GBT visibility data were combined 
and imaged with MIRIAD. The combined images keep 
the native channel width of the VLA data. Image prop¬ 
erties are listed in Table 2. 

We also observed the NH 3 (3,3) transition using the 
VLA with the WIDAR correlator in its D configuration 

10 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 

11 http://gbtidl.nrao.edu/ 
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towards G28.53 MM1 in 2010 May. The calibrators are 
listed in Table 1. The data were calibrated and imaged 
using CASA. Image properties are listed in Table 2. 

2.3. VLA maser observations 

The H 2 0 maser at 22 GHz and the class I CH 3 OH 
maser at 25 GHz were observed using the VLA with the 
WIDAR correlator in 2010 November, toward two posi¬ 
tions, MAI and MA5 which are close to MM1 and MM5. 
Details of observations are listed in Table 1. The data 
were calibrated and imaged using CASA. Image proper¬ 
ties are listed in Table 2. 


3. RESULTS 


3.1. NH 3 emission and temperature 

We detected NH 3 (1,1) and (2,2) emission with both 
VLA and GBT. In the combined images, we found 
NH 3 (1,1) and (2,2) emission associated with from MM1 
to MM8, enabling us to derive temperatures of these 
clumps. The rotation temperature derived from these 
two transitions increases almost linearly with the kinetic 
temperature and differs by <3 K up to 20 K (Walmsley & 
Ungerechts 1983), therefore can be used an equivalence 
of the kinetic temperature. 

We simultaneously fitted the NH 3 (1,1) and (2,2) spec¬ 
tra from the combined images, to derive the best fitted 
linewidth, line intensity, and optical depth at the same 
time. The fitting model includes three gaussians for the 
(1,1) main component and inner satellite components, 
and one gaussian for the (2,2) main component, which 
are all in the form of 


I(v) = I 0 


1 — exp 


-r 0 exp 


/_ 1 ( V-Vp 

( 2\ (T„ 



( 1 ) 


in which Iq is the effective line intensity, tq is the optical 
depth, vq is the central velocity, and a v is the velocity 
dispersion. The central velocity was determined using a 
cross-correlation between the spectra and a model spec¬ 
trum and was then fixed during the fitting. We put con¬ 
straints so that the ratio of optical depths of the main 
and satellite lines of NH 3 (1,1) is 1/0.278 for both satel¬ 
lite lines, all four gaussians have the same velocity dis¬ 
persion, and the velocity separations between the main 
and two inner satellite lines are 7.47 and 7.57 kms' 1 
(Mangum et al. 1992), then we minimized the difference 
between the spectra and the model, using the Levenberg- 
Marquardt algorithm implemented by the lmfit Python 
package 12 . The rotation temperature and NH 3 column 
density were then derived from the best-fitted parame¬ 
ters as in Ho & Townes (1983) and Mangum et al. (1992). 
The temperature map is shown in Figure 2a. The ro¬ 
tation temperature throughout G28.53 is <17 K. Even 
around the dust core ‘core 2’ and the masers, it is not 
showing increased temperature. 

The NH 3 spectra within MM1 present two velocity 
components, at ~85 kms -1 and ~87 kms -1 respectively. 
Similar multiple components have been detected in the 
central regions of a few filamentary IRDCs (Zhang & 
Wang 2011; Henshaw et al. 2013). For these spectra, 
the general fitting procedure above derived overall prop¬ 
erties which blended the two components and derived a 


12 https://github.com/newville/lmfit-py 


temperature that is usually a few K different from the 
true temperature of either component (cf. Figure 2a and 
Figure 3). 

We attempted to fit the two velocity components si¬ 
multaneously and derive a temperature map for each 
component. However, it was unsuccessful thanks to the 
limited signal-to-noise ratios in the spectra. Thus, we 
averaged the spectra within the cores in MM1 defined in 
Section 3.2 to increase the signal-to-noise ratio, and fit¬ 
ted them with a model containing two components each 
of which includes four gaussians as in the single velocity 
component case. The two central velocities were both 
fixed during the fitting. The fitting results are shown in 
Figure 3 and Table 3. Typical temperatures are 13-17 K, 
while the temperature of the 85 kms -1 component is in 
general 2 K higher that of the 87 kms -1 component, sug¬ 
gesting that there are indeed two distinct components. 
We obtained NH 3 column densities of the two velocity 
components, which were used to determine the ratio of 
dust masses in the two components in Section 3.2. 

We also derived the rotation temperature using the 
VLA NH 3 data only. Without GBT data, the interfer¬ 
ometer filters out extended emission and has more con¬ 
tribution from compact structures. However, the change 
of the temperature map is insignificant, as shown in Fig¬ 
ure 2b. The highest temperature is >20 K which is found 
in MM1, while typical temperatures are still 13-17 K as 
in Figure 2a. Note that the two velocity components also 
show up in MM1 in the VLA data. We fitted the aver¬ 
age spectra of MMl-pl with two velocity components 
and showed the result in Figure 4. The temperature of 
the 87 kms -1 component is consistent with that derived 
from the combined data, while the temperature of the 
85 kms -1 component is up to ~20 K. 

In addition, we detected NH 3 (3,3) emission in MM1 
with VLA. Unlike NH 3 (1,1) and (2,2), the (3,3) emission 
is concentrated. We marked the three positions exhibit¬ 
ing strong (3,3) emission in Figure 2b as T1-T3, as well 
as ‘core 2’ which presents weak (5<r) (3,3) emission, with 
their spectra shown in insets. T1 presents two velocity 
components, at ~85 and ~87 kms -1 . The other three 
spectra all center at ~87 kms -1 . We fitted a single gaus¬ 
sian to the spectra of ‘core 2’, T2, and T3, and two gaus¬ 
sians at 85.5 and 87.9 kms -1 to the spectrum of Tl. The 
FWHM linewidth of ‘core 2’ is up to 4.0 kms -1 , while 
those of T2 and T3 are 0.5 and 1.2 kms -1 , respectively. 
The FWHM linewidths of the two components of Tl are 
2.6 and 1.8 kms -1 respectively. It is worth noting that 
T2 and T3 are next to the H 2 0 maser W4, whose posi¬ 
tions are consistent with the possible CO outflow lobes 
surrounding W4 as discussed in Section 3.4. These two 
NH 3 (3,3) emission peaks might suggest outflow heating 
and turbulence injection given the high excitation tem¬ 
perature of NH 3 (3,3) (123.5 K) and the large linewidths 
we detected. This is similar to what Wang et al. (2012) 
found in the IRDC G28.34. 

3.2. Fragmentation in clumps 

The SMA observations of the four clumps revealed 
0.1 pc scale fragmentation. All peaks with fluxes larger 
than 5<r in the SMA 230 GHz dust emission images were 
identified. These ^0.1 pc scale peaks are in agreement 
with the definition of cores (e.g. Zhang et al. 2009). Then 
we fitted 2D gaussian functions to them and obtained 
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their positions and deconvolved sizes. We also applied 
the primary beam correction to the images and obtained 
the corrected fluxes of the cores. The results are listed 
in Table 4. 

We assumed a gas-to-dust mass ratio of 100, and ther¬ 
mal equilibrium between dust and gas, then the mass is 
(Beuther et al. 2005): 


Mr. 


2.0 x 10 -2 F v / D \ 2 / v \ -3-/3 
Jv (T ro t) V1.2 THz/ 


M 0) 

( 2 ) 


in which J^(T rot ) = [exp ( hv/kT mt ) — 1] 1 , and the dust 
emissivity index /? = 1.5 for MM1 and MM8, 1.7 for 
MM5, and 1.0 for MM7 (Rathborne et al. 2010). 

We also reported errors of M core in Table 4, assuming 
that M cole is only dependent on T rot , F v , and D. The 
error in the kinematic distance was assumed to be 10% 
(Reid et al. 2009). The errors in the fluxes and tem¬ 
peratures listed in the table are from gaussian fittings 
without accounting for observational uncertainties. The 
observational uncertainty of the fluxes are mainly due to 
flux calibration as described in Section 2.1, which is typ¬ 
ically 20%, while that of the temperatures depends on 
the signal-to-noise ratio of the NH3 lines (e.g. Li et al. 
2003), which is typically <0.5 K for our spectral data 
with a signal-to-noise ratio of >20. Accumulating all the 
errors listed above, typical errors in M core are 30%, which 
should be interpreted as lower limits given the unknown 
uncertainty in /3, possible discrepancy in the tempera¬ 
tures of dust and NH3, and the uncertainty in gas-to-dust 
mass ratio. 

For MM1, using the compact configuration data, we re¬ 
solved it into six cores at an angular resolution of ~2.3", 
labelled as MMl-pl to MMl-p6 in Figure 5a. We de¬ 
tected CH3OH line emission at 87 kins -1 towards MMl- 
pl, which is the only core that this line is detected among 
all cores we identified in the four clumps. 

The cores in MM1 all exhibit two velocity components 
in NH3 lines. NF13 emission of both components mor¬ 
phologically matches the dust emission. We considered 
the component at ~87 kms -1 , because both the CH3OH 
line emission we detected and the single-dish N 2 H+ and 
H 13 CO + lines (Rathborne et al. 2008) are at ~87 kms -1 . 
After deriving masses using the parameters listed in Ta¬ 
ble 4, we determined the ratios of masses in the two com¬ 
ponents based on the NH 3 column densities in Table 3, 
and scaled the masses of the cores. 

Among the six cores in MM1, MMl-pl is consistent 
with ‘core 2’ reported in Swift (2009) while MMl-p2 is 
consistent with the core reported in Rathborne et al. 
(2008), both using 345 GHz dust continuum emission. 
Rathborne et al. (2008) resolved MMl-p2 into three con¬ 
densations, at an angular resolution of ~2". The three 
condensations have sizes of ~0.06 pc and masses of 9- 
20 Mg. In addition, in the dust emission map of Swift 
(2009), we found weak features at 3-5 ct levels that are 
coincident with the other four cores in MM1. 

For MM5, we identified three cores, which are labelled 
as MM5-pl to MM5-p3 in Figure 5b. MM7 and MM8 
are at the southeast corner of G28.53, in a filamentary 
structure isolated from the main part of the cloud. We 
resolved MM7 to one core, and MM8 to two cores, la¬ 
belled in Figure 5c. The properties of these cores are 


listed in Table 4. 

Aside from dust continuum emission, we detected a 
few spectral lines toward these cores using the SMA. As 
shown in Figure 6, MMl-pl exhibits 12 CO, 13 CO, and 
C 1S 0 2-1 lines, as well as CH 3 OH8(—l,8)-7(0,7) and 
H 2 CO 3(1,2)—2(1,1) lines. MM5-pl exhibits 12 C0 2-1 
and H 2 CO 3(0,3)-2(0,2) lines, but little 13 CO, C 18 0 or 
CH3OH line emission. MM7-pl and MM8-pl exhibit 
only 12 CO line. CH3OH and H 2 CO lines are tracers of 
dense gas given their large critical densities (10 5 6 cm -3 ), 
while 12 CO line could trace outflows. The relationship 
between these lines and star formation in the cores will 
be discussed in Section 4.3. 

3.3. Protostellar outflows 

The SMA observations covered the 12 C0 2-1 line, 
which can trace molecular outflows driven by protostars. 
We detected CO outflows in MM1 and MM5. MM7 and 
MM8 also present CO emission, but do not show any 
signature of outflows. 

The CO channel maps in Figure 7 reveal two parallel 
bipolar outflows, associated with MMl-pl and MMl-p2. 
The CO emission associated with MMl-pl presents two 
blue shifted components at 81 and 83 kms -1 with re¬ 
spect to the system velocity of 87 kms -1 , on opposite 
sides of MMl-pl. At more redshifted velocities between 
90-100 kms -1 , CO emission is also found around MMl- 
pl. It is consistent with a bipolar outflow nearly parallel 
to the plane of the sky, with the two components at 81 
and 83 kms -1 tracing its near side while the emission at 
more redshifted velocities tracing its far side (e.g. Fig¬ 
ure 8 of Wu et al. 2009). Similarly, the CO emission 
associated with MMl-p2 reveals two components at 90- 
92 kms -1 and 92-94 kms -1 on two opposite sides of 
the core, which might trace the red-shifted component 
of a bipolar outflow in the plane of the sky. We found a 
CO emission gap at 85-88 kms -1 , indicating the SMA 
filtered out the extended CO emission around the sys¬ 
tematic velocity. 

A H 2 0 maser detected in Wang et al. (2006) (Wl, 
Section 3.4) is in the projected path of the CO outflow 
associated with MMl-p2. Its velocity is 85 kms -1 , while 
the CO velocity here is 94-96 kms -1 which is 10 kms -1 
apart, therefore it is unlikely excited by the outflow. We 
also found compact H 2 CO line emission, which has a 
critical density of ~10 5 cm -3 (Guzman et al. 2011), at the 
same position of Wl but at the velocity of 89 kms -1 . It 
could be a dense core which is not detected by the SMA 
dust observations because of insufficient sensitivity at the 
edge of the primary beam. 

We also found a bipolar outflow in MM5, shown in the 
CO channel maps in Figure 8. The outflow is oriented 
in a north-south direction, with a blue shifted compo¬ 
nent at 79-83 kms -1 and a red shifted component at 
90-96 kms -1 . The central velocity, ~87 kms -1 , is con¬ 
sistent with the velocities of the cores. However, the 
three cores in MM5 are closely packed, thus we can not 
identify which one is driving this outflow. 

Assuming LTE and optically thin 12 CO emission in the 
outflow wings, we derived outflow properties as in Wang 
et al. (2011), listed in Table 5. The relative abundance 
[CO]/[H 2 ] was assumed to be 10 -4 (Blake et al. 1987). 
We did not correct for inclination angles of the outflows. 
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The dynamic ages of ^10 4 years suggest very early 
evolutionary phases of these cores (cf. Zhang et al. 2001, 
2005; Beuther et al. 2002). Among the three outflows, 
the two in MM1 are more massive and more energetic 
than the one in MM5. All three outflows are comparable 
with the high-mass outflows in other IRDCs (e.g. Wang 
et al. 2011, 2014). 


3.4. H 2 O and CH 3 OH masers 

The VLA observations of class I CH3OH maser at 
25 GHz did not detect any sources. The 22 GHz VLA ob¬ 
servations revealed three H2O masers around MM1, and 
none in MM5. Including the H2O maser found by Wang 
et al. (2006), there are 4 H2O masers around MM1, la¬ 
belled as W1-W4 in Figures 2 and 5a, among which W1 
and W4 are within MM1. Table 6 summarizes properties 
of the masers. 

H2O masers in star forming regions are believed to 
be excited in shocked ambient gas (Elitzur et al. 1989), 
therefore could trace protostellar outflows. Surveys have 
found H2O masers to be associated with protostars with 
a wide range of luminosities (e.g. Furuya et al. 2001; 
Szymczak et al. 2005). H 2 0 masers have been used as 
star formation indicators in IRDCs (Wang et al. 2006). 
As discussed above, W1 is not likely associated with the 
outflow driven by MMl-p2, but is coincident with a dense 
core revealed by H2CO. It is also associated with an in¬ 
frared source as shown in Figure 1. W2 is to the north 
of MM1, far away from the field of the SMA observa¬ 
tions, so we could not relate it to any dense cores or 
outflows. It is spatially coincident with MM9 defined 
by Rathborne et al. (2006). We also found NH 3 emis¬ 
sion as well as a 24 /ini infrared source associated with 
it in Figures 1 and 2. W3 is also beyond the scope of 
the SMA observations, spatially coincident with MM4 
defined by Rathborne et al. (2006). A dust core found 
by Swift (2009) is likely coincident with it, although the 
association is uncertain because of a lack of measure¬ 
ment of the core velocity. W4 is within the SMA field, 
and is coincident with weak dust emission as shown in 
Figure 5a. In Figure 7, there is a strong CO emission 
component to the south of W4 at 90-100 kms -1 , while 
only weak CO emission shows up symmetrically to the 
north of it at 81 kms -1 , resulting in a central velocity 
of ~85-86 kms -1 if there was an outflow, which is con¬ 
sistent with the velocity of W4. Therefore, W1 and W4 
might trace two protostars within MM1 that are missed 
by the SMA continuum observations. 


4. DISCUSSIONS 
4.1. Gas fragmentation 

The SMA observations have resolved the four clumps 
into 12 cores in total in dust emission. We compare our 
result with what Jeans fragmentation predicts, in which 
an initially homogeneous piece of gas will fragment into 
smaller pieces with typical size of the Jeans length and 
typical mass of the Jeans mass, given a characteristic 
velocity and a particle number density. In our case, the 
four clumps fragment into cores, thus we can calculate 
the Jeans masses of the clumps and compare them with 
the core masses. 


The Jeans length is defined as 


Aj = v 


- 


1/2 


GpJ 

°' n PC (o.l kms -1 ) ( 


10 4 


cm 


-3 


-1/2 


( 3 ) 


where v is either the sound speed c s or the velocity dis¬ 
persion a v , and n is the particle density of the clump. 
We used the mean NH 3 temperature of the clump to cal¬ 
culate c s . With the parameters reported in Rathborne 
et al. (2006, 2010), we derived the density of each clump 
as n = mass/(4/3-7T-radius 3 ) x l/(pmn), where the mean 
molecular mass number is p = 2.33. Note that for MM5, 
we used the 1.2 mm flux in Rathborne et al. (2006), the 
dust emissivity index in Rathborne et al. (2010), and the 
mean NH 3 rotation temperature to calculate the clump 
mass using Equation 2, instead of directly quoting the 
mass in Rathborne et al. (2010), because of the possi¬ 
bly incomplete spectral energy distribution (SED) fit¬ 
ting (see Section 4.3.3). The resulting clump mass is 
4 times larger. For consistency, we calculated masses of 
the other three clumps in the same way, and their masses 
are within a factor of 50% with respect to the results of 
Rathborne et al. (2010). 

The Jeans mass is then the mass enclosed by a sphere 
with a radius of Aj/2: 


Mj 


47T p / Aj \ 
3 \2j 


= 0.43 M 0 f-— — j-) 3 ( U ) 

V0.1 kms V VlO 4 cm 3 / 


-1/2 


( 4 ) 


The Jeans masses of the clumps are listed in Table 4. 

If the sound speed c s is used to calculate the Jeans 
parameters, which indicates that thermal pressure dom¬ 
inates during fragmentation, the Jeans masses are much 
smaller than the core masses. When the observed ve¬ 
locity dispersion a v is applied in Equation 4, the Jeans 
masses are then consistent with the observed core masses. 
Therefore, turbulent pressure is dominant over thermal 
pressure during 1—>-0.1 pc scales fragmentation, which 
is consistent with fragmentation studies of other IRDCs 
(Zhang et al. 2009; Zhang & Wang 2011; Wang et al. 
2011, 2014). 


4.2. Gravitational collapse of cores 

With the velocity dispersion a v and the deconvolved 
radius R based on Table 4, we analyze virial status of 
the cores. The gas is in a virial equilibrium if the kinetic 
energy equals half of the gravitational energy. If the core 
mass is larger than its virial mass, it is gravitationally 
bound and will collapse. Note that in this framework, 
rotation of cores, magnetic field, or external pressure are 
all ignored, therefore the virial analysis is only robust in 
the simplest scenario. 

The virial mass depends on density profile of the core 
(e.g. MacLaren et al. 1988). Assuming the radial density 
profile is p[r) ~ r ~ 2 as measured in IRDC cores (Wang 
et al. 2011), then the virial mass is 


Afvirial — —^— — 69.8 Mg 


( _ ° v ) 2 ( R ^ 

VI kms 1 / \ 0.1 pc/ 
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If instead the density does not depend on the radius, 
p[r) = const, then will be 1.7 times larger. The 

viral parameter is defined as a = M v i r i a i/M core . The 
result is listed in Table 4. 

All of the 12 cores have a virial mass smaller than the 
core mass (a < 1). The virial parameters are between 
0.09 and 0.7 in general. MMl-pl has the smallest a of 
0.09. Even assuming a uniform density profile, all but 
MMl-p3, MM7-pl, and MM8-p2 have a < 1. 

Therefore, in the simplest scenario where only grav¬ 
itational and kinetic energies are considered, all cores 
are gravitationally bound and tend to collapse if p{r) ~ 
r~ 2 . The cores with significant star formation, MM1- 
pl/p2 and MM5-pl (see Section 4.3), are strongly self- 
gravitating (a < 0.2). Even in the relative quiescent 
clumps such as MM7 and MM8, the cores are still grav¬ 
itationally bound. 

We compared the virial parameters of the cores in 
G28.53 with those of hot cores in the NH 3 sample in our 
recent work by Lu et al. (2014) in Figure 9. The three hot 
core sources are IRAS 18089-1732, IRAS 18360-0537, 
and IRAS 18414—0339. Note that in Lu et al. (2014), the 
core masses are determined from the NH 3 column den¬ 
sities, although for the ones that dust continuum data 
are available the masses based on the two datasets are 
fairly consistent. The rotation temperatures of the hot 
cores measured from the first two inversion transitions 
of NH 3 should be treated as the lower limits, while the 
real temperatures at central part of the hot cores could 
be >300 K. 

Assuming a uniform radial density profile, the virial 
parameters of the hot cores are all smaller than 1, sim¬ 
ilar to what we find in G28.53. The gas temperatures 
of the hot cores are usually >30 K, thus their thermal 
linewidths are larger than those of G28.53. However as 
shown in Figure 9 the non-thermal linewidths are always 
dominant over the thermal linewidths, thus the increase 
of linewidths by a factor of 2-3 from G28.53 cores to hot 
cores is mainly due to the non-thermal motions. There¬ 
fore, the fact that the virial parameters are consistent 
from G28.53 cores to hot cores indicates that turbulence 
is essential in supporting the cores against gravity con¬ 
stantly in the evolution of dense cores. Given that the hot 
cores are more massive than the IRDC cores (typically a 
few hundreds of M 0 ), to maintain a similar virial param¬ 
eter, turbulence must be enhanced. The feedback from 
protostars in hot cores, including outflows and heating, 
could increase the gas temperatures, and more impor¬ 
tantly, inject turbulent energy into gas, which maintains 
a steady accretion of dense gas to form high-mass stars. 

As one of the most important energy sources against 
gravity other than kinetic energy, magnetic field is not 
considered since we do not have measurements in this 
cloud. If the typical magnetic field strength B = 1 mG 
in high-mass star formation regions (Cortes & Crutcher 
2006; Girart et al. 2013; Qiu et al. 2013, 2014; Frau et al. 
2014; Zhang et al. 2014) is added into the virial relation, 
and assuming a uniform radial density profile, then 

(3a 2 + iai)M-^^ = 0, (6) 

in which a a = B/y/4np is the Alfvenic velocity. The 


total virial mass M is 

M = M virial + M b = 
where M v i r j a i = 5afR/G and 
M b = 5a 2 A R/6G = 


5(cr 2 + cri/6)i? 


G 


(M, 


mag 


\ A/core 

is the magnetic virial mass. Here 

R 2 B 


M, 


mag 


M mag = 


V/18G/5 


(7) 

( 8 ) 

(9) 


is the traditionally defined critical mass for a spherical 
cloud of uniform density with a uniform magnetic field 
(Strittmatter 1966). Only when M mag = M cole does the 
magnetic virial mass Mb we defined equal M mag . 

In this case, the virial parameters of MMl-pl/p2 and 
MM5-pl are still smaller than 1, while all the other cores 
have a > 1. This might explain why signatures of star 
formation are only detected in these three cores (see Sec¬ 
tion 4.3). Therefore, magnetic field might also be an im¬ 
portant factor in determining the dynamics of the cores. 


4.3. Massive star formation in MM1 and MM5 
4.3.1. An outflow-disk system in MMl-pl 

In several previous studies (Rathborne et al. 2008, 
2010; Swift 2009; Sanhueza et al. 2012), MM1 was clas¬ 
sified as ‘quiescent’ with no signs of high-mass star for¬ 
mation, given the lack of molecular lines and absence 
of associated infrared sources. Our observation suggests 
that MM1 is not quiescent, but is already forming stars. 
In MM1, the two most massive cores are associated with 
outflows. The two H2O masers W1 and W4 could also 
trace star formation. In total there are at least four star 
forming sites. 

With a mass of 56.0 M@, MMl-pl is the most massive 
core we resolved in G28.53. It remains unresolved in 
the 345 GHz SMA observations at an angular resolution 
of 2" x 1.1" (Swift 2009), and only one 12 CO outflow 
was detected, both of which suggest there might be a 
single protostellar object embedded in MMl-pl, which 
will form a star or a multiple stellar system. The large 
core mass and the energetic outflow suggest high-mass 
star formation. The small virial parameter of 0.09-0.2 
suggests that it is strongly self-gravitating. 

The SMA CH 3 OH line emission at ~87 kms” 1 in 
MMl-pl presents a flatten morphology perpendicular to 
the CO outflow, with a velocity gradient of ^3 kms -1 
across ^0.1 pc, shown in Figure 10a. With a critical den¬ 
sity of ~10 5 cm -3 , the CH 3 OH line emission is tracing 
a dense gas component, likely a rotating envelope sur¬ 
rounding a protostellar disk. These results suggest that 
there might be an outflow-disk system within MMl-pl. 


4.3.2. Star formation activity associated with MMl-p2 and 
cluster formation in MM1 

As Rathborne et al. (2008) pointed out, MMl-p2 frag¬ 
ments into three condensations among which the most 
massive one might form high-mass stars, based on the 
broad line wings in singled-dish spectra of the outflow 
tracer SiO. With the SMA 230 GHz observations we de¬ 
tect the outflow directly, although it is uncertain which 
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one of the three condensations in Rathborne et al. (2008) 
is driving this outflow. 

The outflows associated with MMl-pl/p2 are parallel 
to each other. Similar nearly parallel outflows have been 
found in other IRDCs (e.g. Wang et al. 2011, 2014). A 
possible interpretation is that the cores are fragmented 
from an initially rotating clump, then their rotation axes 
will tend to be parallel to the initial rotation axis of the 
clump given conservations of angular momentum, which 
leads to parallel disks and finally parallel outflows. As 
shown in Figure 11, the velocity field of the 87 kms ' 1 
component of GBT NH 3 (2,2) of MM1 presents a gra¬ 
dient at a position angle of 40-50° north to east, per¬ 
pendicular to the outflows. If this velocity gradient is 
interpreted as rotation, then the MM1 clump is likely 
spinning slowly, with an axis parallel to the outflow ori¬ 
entations, which supports the rotating clump fragmenta¬ 
tion scenario. However, note that such velocity gradient 
can also be reproduced with random motions (Burkert 
& Bodenheimer 2000). In addition, magnetic field in the 
clumps might also be important in aligning the outflows 
(e.g. Zhang et al. 2014), by regulating the 0.1 pc scale 
magnetic fields in the cores thus affecting the orientation 
of the accretion disks. 

The two H 2 O masers in MM1 could pinpoint another 
two star forming sites. Both are associated with weak 
(3cr) 345 GHz dust emission in the SMA map of Swift 
(2009). Between them, W1 is more interesting because 
of its association with H 2 CO line emission and an in¬ 
frared source. Three more cores were detected at levels 
of > 6 cr in the southern part of MM1 in the SMA 345 GHz 
dust emission map of Swift (2009). Therefore, with at 
least eleven cores, among which four are forming stars, 
MM1 will probably form a star cluster. MMl-pl and the 
most massive condensation in MMl-p2 are two promising 
high-mass star formation sites. 

In Figure 1, MM1 is the most infrared extinct part 
in G28.53, and it remains dark up to 70 pm (Molinari 
et al. 2010), which leads to a low luminosity (e.g. 200 Lq 
in Rathborne et al. 2010). In contrast, MM5 starts to 
be bright in 24 pm and has a larger luminosity (Rath¬ 
borne et al. 2010). The discrepancy between the more 
active star formation and the smaller luminosity of MM1 
could be due to heavy extinction: if the 85 kms -1 com¬ 
ponent was at the front then it could extinct infrared 
emission produced by star formation in the 87 kms ' 1 
component behind it. Even if the 87 kms -1 component is 
not shielded, the self-extinction ( Ay up to 50-100 given 
a column density of >10 23 cm -2 , Giiver & Ozel 2009) 
could keep it infrared dark. This is similar to the IRDC 
G30.88+0.03 in Zhang & Wang (2011), which has a small 
luminosity of 460 Lq but presents H 2 O masers and mas¬ 
sive cores, and there are also two velocity components 
along the line of sight. Therefore, infrared emission only 
is not sufficient in assessing star formation activity, es¬ 
pecially in the presence of multiple velocity components, 
while spectral lines like CH 3 OH can be used to resolve 
multiple components and trace deeply embedded star for¬ 
mation. 

4.3.3. Early phase star formation in MM5 

Previous studies (Rathborne et al. 2010; Sanhueza 
et al. 2012) classified MM5 as ‘intermediate’, based on 


infrared emission and chemistry, e.g. it contains a 24 pm 
source but not a green fuzzy at 4.5 pm. Our results 
support this conclusion, with MM5 in an ‘intermediate’ 
evolutionary phase between the ‘active’ MM1 and the 
‘quiescent’ MM7/8. 

A bipolar outflow in the north-south direction is found 
in MM5, while the angular resolution is not good enough 
to tell which of the three cores it originates from. All 
three cores are self-gravitating, with masses between 24 
and 38 M 0 , thus are equally possible to harbor proto¬ 
stars. The outflow is less energetic and slightly younger 
than the two outflows in MM1, suggesting that the star 
formation in MM5 is less active and possibly at an earlier 
phase than in MM1. 

With a critical density of >10 5 cm -3 (Guzman et al. 
2011), H 2 CO usually traces dense gas. As shown in Fig¬ 
ure 10b, the SMA detected H 2 CO line emission associ¬ 
ated with MM5-pl/p2, which is elongated in the east- 
west direction and presents a velocity gradient across 
MM5-pl and MM5-p2, perpendicular to the outflow. 
Therefore, the H 2 CO emission is likely tracing a dense 
envelope that encloses the two cores, and perhaps spins 
slowly, which might feed material into disks around em¬ 
bedded protostars. The two cores, MM5-pl/p2, are 
hence the potential driving sources of the outflow. 

The NH 3 rotation temperatures in MM5, using either 
combined data or VLA data only, are generally <15 K. 
We did not find any signature of heating associated with 
outflows or hot cores. The protostars in the three cores, if 
there are any, must be at an early stage when accretion 
or stellar luminosity is not strong enough to heat the 
ambient gas. Note that Rathborne et al. (2010) derived 
a dust temperature of 30 K and a luminosity of >10 3 L 0 
for MM5, from SED modeling. The robustness of this 
result suffers from the fact that they used fluxes in only 
four bands (24 pm, 450 pm, 850 pm, 1.2 mm), making 
the modeling less constrained at around 100 pm which is 
supposed to be the peak of the SED. Recently Schneider 
et al. (2014) fitted the SED of G28.53 using Herschel 150- 
500 pm data and derived a dust temperature of 16-18 K 
for MM5, consistent with our NH 3 rotation temperature, 
which is expected since at high densities (> 10 4 ’ 5 cm -3 ) 
the gas and dust temperatures are usually well coupled 
(Goldsmith 2001). The results of MM1 and MM7 which 
include data from more bands in Rathborne et al. (2010) 
are in general consistent with ours. 

With a total mass of ~400 M 0 (Rathborne et al. 2006, 
2010), MM5 has sufficient dense gas to form ~80 M 0 
stars assuming a star formation efficiency of 20%. The 
24 pm point source, the outflow, the spinning H 2 CO en¬ 
velope, and the self-gravitating cores all suggest the ca¬ 
pability of liigh-mass star formation, while the absence of 
hot core tracers such as CH 3 OH and the low temperature 
indicate little stellar feedback. 

4.4. Potential of star formation of MM7 and MM8 

The three cores in MM7 /8 are marginally gravitation¬ 
ally bound, and do not show any star formation signa¬ 
ture such as outflow, maser, infrared source or increased 
temperature. As previous studies suggested (Rathborne 
et al. 2010; Sanhueza et al. 2012), these two clumps are 
likely ‘quiescent’, with no current star formation. 

However, given that the cores are self-gravitating, 
MM7 /8 still have the potential of collapsing and eventu- 
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ally forming stars. Their masses are ~300 and ~240 M 0 
respectively using the fluxes and /3 reported in Rathborne 
et al. (2010), which are comparable to the mass of MM5. 
With mean densities of 1 x10 4 -2x10 4 cm -3 , the free fall 
time scale is 0.2-0.4 Myr. If the cores are indeed collaps¬ 
ing even after taking magnetic field into account, then 
star formation will start in such a time scale. 

It is also worth noting that MM7 /8 are embedded in a 
filamentary structure, which presents a velocity gradient 
of 1 kms ' 1 along the major axis, which is consistent 
with a gas stream falling into the main part of G28.53, 
evidenced by the GBT NH 3 (2,2) velocities in Figure 11. 
In fact, G28.53 is the densest and coldest part of a large 
filament which extends ~60 pc (Wang et al. submitted), 
in which MM7/8 connect the main part of G28.53 and 
the large filament. Assuming a proper motion velocity of 
~l-2 kms -1 and a projected distance of 2-3 pc, MM7/8 
will collide with MM1 in 1-3 Myr. For reference, the free 
fall time scale of the main part of G28.53 (MM1/2/3/4/6, 
~3.5xl0 3 Mq within 3 pc; Rathborne et al. 2010) is 
~1.4 Myr. Therefore MM7/8 will probably first form 
stars, then dynamically interact with the other clumps at 
the center of G28.53 due to their gravitational pull. This 
picture is similar to what Liu et al. (2012) proposed for 
the high-mass cluster formation in G10.6—0.4, in which 
fragmentation occurs in filaments on free fall time scales, 
while the filaments themselves fall into a central dense 
region where OB clusters are forming through a global 
contraction. 

5. CONCLUSIONS 

The multi-frequency observations of the IRDC G28.53 
have revealed twelve 0.1 pc scale cores in four massive 
clumps. We obtain gas temperatures and kinematics 
from VLA and GBT NH 3 lines, and derive core masses 
based on the SMA dust continuum emission. We also 
assess star formation associated with these cores using 
VLA masers and the SMA spectral lines. Our conclu¬ 
sions are: 

1. When 1 pc scale clumps in IRDCs fragment into 
0.1 pc scale cores, turbulent pressure is more im¬ 
portant than thermal pressure in supporting the 
gas. Therefore, it is possible to form super- 
thermal-Jeans-mass cores that are massive enough 
to harbor high-mass stars. 

2. All of the cores we find in IRDC G28.53 have a 
virial parameter a (M v i r i a i/M cole ) smaller than 1 
assuming a radial density profile of p(r ) ~ r^ 2 , and 
nine out of twelve cores still have a < 1 even assum¬ 
ing a constant radial density, suggesting that they 
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are all gravitationally bound, if magnetic field, ro¬ 
tation, or external pressure are not considered. The 
virial parameters of the three star forming cores in 
MM1 and MM5 are -C 1, therefore these cores are 
strongly self-gravitating. Even if magnetic field is 
taken into account, these three cores are still grav¬ 
itationally bound, consistent with star formation 
activities associated with them. 

3. The virial parameters of the cores in G28.53 and 
in several hot core sources are all ^ 0 . 2 - 1 , even 
though gas temperatures of hot cores are much 
higher, suggesting turbulent motions are more im¬ 
portant than thermal motions in supporting mas¬ 
sive cores against gravity as they collapse. The con¬ 
sistent virial parameters from the lower-mass cores 
in IRDCs to the more massive hot cores suggest 
that turbulence is enhanced as cores are accreting, 
which is essential in maintaining a steady accumu¬ 
lation of dense gas thus forming high-mass stars. 

4. The most massive core in MM1 is found at the 
center of this clump. The energetic outflow traced 
by CO and a disk-like structure traced by CH 3 OH 
suggest that there might be an outflow-disk system 
in this core. With a mass of 56 M Q , it will form a 
high-mass star or multiple stellar system. 

5. In addition, MM1 harbors another three possible 
star forming sites, as well as seven cores that are 
likely quiescent but also massive. Therefore MM1 
might be a progenitor of high-mass star cluster, 
rather than a quiescent clump suggested by previ¬ 
ous studies. 

6 . In contrast, MM5 clump is intermediate in the evo¬ 
lutionary stage and MM7/8 clumps are quiescent. 
However they have the potential to continue frag¬ 
ment and collapse until star formation starts as in 
MM1. MM7/8 might be falling into the main part 
of G28.53 due to its gravitational pull while frag¬ 
mentation occurs in themselves, which provides a 
way to form star clusters. 
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TABLE 1 
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Jonathan Swift 

2009A-H004 

H 2 co 

2009 06 07 

MMl 

3C273 

Callisto 

1751+096 

SMA compact 

Ke Wang 

2010A-S053 

12 co, 13 co, c 18 o, h 2 co 

2010 06 20 

MM5, MM7, MM 8 

3C454.3 

Titan, Neptune 

1911-201 

VLA D 

Ke Wang 

AW763 

NH 3 (1,1), (2,2) 

2010 01 04 

A3, A4, A5 

3C373 

3C286 

1851+005 

VLA D 

Ke Wang 

AW763 

NH 3 (1,1), (2,2) 

2010 01 08 

Al, A2 

3C454.3 

3C48 

1851+005 

VLA D 

Ke Wang 

AW776 

NH 3 (3,3) 

2010 05 09 

MMl 

3C454.3 

3C48 

1851+005 

GBT 

Ke Wang 

AGBT10A-067 

NH 3 (1,1), (2,2) 

2010 02 28 

8 ' x7.5' region 

3C353 

3C48 

1852+0055(ptg) 

VLA C 

Ke Wang 

10B-182 

H 2 CO, CH 3 OH 

2010 11 24 

MAI, MA5 

3C454.3 

3C48 

1851+005 


a Coordinates of pointing centers: MM1: (18:44:18.0, -3:59:23.0); MM5: (18:44:17.0, -4:02:04.0); MM7: (18:44:23.7, -4:02:09.0); MM8: (18:44:22.5, -4:01:47.0); Al: (18:44:24.0, 
-4:02:12.0); A2: (18:44:17.2, -4:02:00.0); A3: (18:44:16.0, -4:00:42.0); A4: (18:44:18.0, -3:59:35.0); A5: (18:44:18.2, -3:58:40.0); MAI: (18:44:18.26, -3:59:45.50); MA5: 
(18:44:18.30, -4:01:55.00). 
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TABLE 2 
Image properties. 


Image 

Bmaj 

(") 

Bmin 

(") 

Bpa 

(°) 

Band/Channel Width 

RMS 

(mjy beam -1 ) 

SMA MM1 cont. 

3.2 

1.7 

73 

4 GHz 

1.1 

SMA MM5 cont. 

3.2 

2.8 

37 

8 GHz 

0.9 

SMA MM7/8 cont. 

3.1 

2.7 

53 

8 GHz 

0.8 

SMA MM1 spec. 

3.1 

1.6 

72 

1.1 kms -1 

60 

SMA MM5 spec. 

3.2 

2.9 

38 

1.1 kms -1 

50 

SMA MM7/8 spec. 

3.2 

2.9 

39 

1.1 kms -1 

50 

VLA+GBT NH 3 (1,1), (2,2) 

3.9 

2.7 

20 

0.3 kms -1 

3.0 

VLA NH 3 (1,1), (2,2) 

3.8 

2.7 

20 

0.3 kms -1 

3.2 

VLA NH 3 (3,3) 

3.4 

2.6 

-161 

0.2 kms -1 

4.0 

VLA H2O maser 

1.8 

1.0 

34 

0.8 kms -1 

1.1 


TABLE 3 

Fitting results of two velocity components of NH3 spectra in MM1. 


Core ID 

Vi sr a 

(kms -1 ) 

(mJy beam -1 ) 

FWHM 

(kms -1 ) 

T (l,l,m) 

Trot 

(K) 

%h 3 

(10 16 cm -2 ) 

MMl-pl 

85.9 

87.7 

35.0 

25.0 

1.08i0.05 

1.34i0.08 

5.4i0.5 

6.4i0.7 

16.9i0.9 

14.9i0.9 

l.OiO.l 

1.3i0.1 

MMl-p2 

85.2 

88.0 

19.0±1.3 

25.2i0.9 

1.30i0.08 

1.33i0.04 

3.7i0.5 

10.Oil.1 

15.1i0.9 

13.4i0.4 

0.8i0.1 

1.9i0.2 

MMl-p3 

85.1 

87.8 

23.3il.l 

21.8il.3 

1.65i0.07 

1.31i0.07 

6.4i0.6 

5.0i0.6 

13.7i0.5 

14.4i0.7 

1.5i0.1 

l.OiO.l 

MMl-p4 

85.7 

87.8 

28.1il.5 

18.5il.3 

1.24i0.05 

1.24i0.07 

6.6i0.7 

9.4il,8 

15.5i0.6 

13.6i0.8 

1.3i0.1 

1.7i0.3 

MMl-p5 

85.1 

87.9 

17.6il.O 

22.6i0.9 

1.63i0.08 

1.21i0.04 

5.0i0.5 

10.8il.4 

14.4i0.7 

13.6i0.5 

1.2i0.1 

1.9i0.2 

MMl-p6 

85.8 

87.8 

13.0 

19.7il.2 

1.74i0.18 

1.26i0.06 

3.4i0.6 

12.7i2.8 

15.3il.5 

13.1i0.8 

0.9i0.1 

2.2i0.4 


a Velocities of the components were all fixed in the fitting. 

k Intensities of the two components in MMl-pl and one component in MMl-p6 were fixed to make sure the fitting converged. 
















TABLE 4 
Core properties. 


Core ID 

RA 

(J2000) 

Dec 

(J2000) 

(kms x ) 

Maj.xMin. b 

("x") 

PA b 

(°) 

Flux c 

(mjy) 

Trot 

(K) 

FWHM 

(kms -1 ) 

Mj d 

(M@) 

M viria i e 

(M 0 ) 

^virial, B = lmC 

(M 0 ) 

, f M core g 
(M 0 ) 

MMl-pl 

18:44:18.04 

-03:59:22.81 

87.7 

2.2x1.3 

103 

84.04=3.3 

14.94=0.9 

1.34 

1.9 

35.8 

5.7 

9.5 

9.9 

56.04=16.5 

MMl-p2 

18:44:17.76 

-03:59:34.31 

88.0 

5.8x2.2 

147 

46.84=3.6 

13.44=0.4 

1.33 

1.9 

35.8 

11.7 

19.6 

29.5 

46.24=13.7 

MMl-p3 

18:44:17.29 

-03:59:22.74 

87.8 

2 .9x1.1 

22 

15.04=2.1 

14.44=0.7 

1.31 

1.9 

35.8 

5.8 

9.6 

13.5 

7.34=2.4 

MMl-p4 

18:44:18.02 

-03:59:18.96 

87.8 

4.0x1.4 

124 

17.6±0.4 

13.64:0.8 

1.24 

1.9 

35.8 

6.9 

11.5 

18.1 

13.34=3.9 

MMl-p5 

18:44:17.60 

-03:59:29.85 

87.9 

3.8x2.7 

40 

21.04=2.1 

13.64=0.5 

1.21 

1.9 

35.8 

8.9 

14.9 

32.1 

17.24=5.2 

MMl-p6 

18:44:18.06 

-03:59:36.24 

87.8 

3.4x2.0 

89 

16.34=0.2 

13.14=0.8 

1.26 

1.9 

35.8 

7.8 

13.0 

20.9 

16.44=4.8 

MM5-pl 

18:44:17.30 

-04:02:04.88 

87.2 

7.3x0.8 

177 

20.74=2.7 

13.64:0.4 

1.06 

1.4 

21.3 

5.4 

9.0 

11.5 

38.24=12.0 

MM5-p2 

18:44:17.01 

-04:02:01.54 

86.9 

4.4x2.0 

32 

13.54=1.2 

13.84=0.4 

1.30 

1.4 

21.3 

9.4 

15.7 

24.6 

24.34=7.3 

MM5-p3 

18:44:17.16 

-04:02:10.01 

87.0 

5.8x2.8 

40 

14.64=0.1 

12.54=0,4 

1.12 

1.4 

21.3 

9.7 

16.2 

40.6 

30.44=8.7 

MM7-pl 

18:44:24.06 

-04:02:12.60 

88.2 

7.9x2.7 

144 

19.64=4.2 

13.74=0.4 

1.00 

6.0 

58.5 

9.4 

15.6 

130.5 

11.14=4.0 

MM8-pl 

18:44:22.18 

-04:01:43.30 

88.1 

<8.0x2.7 


18.54=1.4 

13.54=0.5 

1.05 

2.9 

37.9 

10.2 

17.0 

69.9 

24.74=7.3 

MM8-p2 

18:44:21.83 

-04:01:38.38 

88.0 

<5.4x2.9 


14.54=0.8 

14.64=0.7 

1.23 

2.9 

37.9 

11.5 

19.1 

58.9 

17.34=5.1 


a Visr of the cores were determined from cross-correlations between a model and the NH 3 (1,1) spectra. For the cores in MM1 that exhibit two velocity components, the ones at 
~87 kms -1 are listed (cf. Table 3). 

k Major and minor FWHMs and position angles of the cores were deconvolved from beam. For MM 8 -pl and MM8-p2 that are not resolved, upper limits of major and minor FWHMs 
are given. 

c Fluxes are corrected for primary-beam response. Errors listed here are from 2D gaussian fittings, not including 20% uncertainty in flux calibration of the SMA data. 

^ The two columns list thermal Jeans masses and turbulent Jeans masses of the clumps, respectively. 

e The two columns list virial masses assuming radial density profile of p ~ r ~ 2 and constant radial density profile respectively. 

^ The total virial mass assuming constant radial density profile and magnetic field strength B —1 mG. 

® For the cores in MM1, the masses have been scaled to those of the 87 kms -1 component, based on the ratios of NH 3 column densities of the two components listed in Table 3. 
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TABLE 5 

CO OUTFLOW PROPERTIES. 


Parameter 

MMl-pl a 

Blue Red 

MMl-p2 b 

Blue Red 

MM5 

Blue Red 

Velocity range (kms -1 ) 

[78,87] 

[78,87] 

[88,96] 

[88,96] 

[ 77 , 88 ] 

[89,99] 

Excitation temperature (K) 

14.9 

14.9 

13.4 

13.4 

13.6 

13.6 

Total mass (Mq) 

0.84 

0.63 

0.71 

1.08 

0.16 

0.24 

Momentum (Mq kms -1 ) 

3.76 

2.66 

3.57 

6.06 

0.58 

1.66 

Energy (Mq km 2 s -2 ) 

10.44 

6.88 

9.81 

18.99 

1.54 

6.55 

Projected lobe length (pc) 

0.52 

0.92 

0.60 

0.71 

0.31 

0.68 

Dynamic age (10 4 yr) 

5.44 

9.52 

7.46 

8.76 

3.10 

5.50 

Outflow rate (10 -5 Mgyr -1 ) 

1.54 

0.66 

0.96 

1.23 

0.53 

0.43 


a Red and blue lobes are both blue-shifted with respect to the core, 
k Red and blue lobes are both red-shifted with respect to the core. 


TABLE 6 

H2O MASER PROPERTIES. 


Maser ID 

RA 

Dec 

Peak flux 

(Jy) 

Peak velocity 
(kms -1 ) 

Velocity range 
(kms -1 ) 

W1 

18:44:16.76 

-03:59:17.00 

6.47 

85.0 

84.4-85.6 

W2 

18:44:17.23 

-03:58:22.50 

1.03 

85.3 

83.4-86.7 

W3 

18:44:18.23 

-04:00:01.50 

0.24, 0.59, 0.18 

83.4, 90.1, 97.7 

81.7-99.4 

W4 

18:44:19.63 

-03:59:32.00 

0.74 

85.9 

84.2-87.6 



Fig. 1.— Spitzer and Herschel composite image overlaid with 1.2 mm continuum emission of G28.53 taken from Rathborne et al. (2006). 
The colors represent emission at 8.0 pm (blue), 24 pm (green), and 70 pm (red) (Churchwell et al. 2009; Carey et al. 2009; Molinari et al. 
2010). The contours levels are 30 (3a), 60, 90, 120, 180, 240 mjy per beam. The four clumps, MM1, MM5, MM7, and MM8 are marked. 
The green circles mark MMl-pl, MM5-1, MM7-pl, and MM8-pl, the most massive cores in each clump. The four crosses mark the H 2 O 
masers. 
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Fig. 2.— NH 3 rotation temperature of G28.53. (a) Temperature map using the combined GBT and VLA data, overlaid with contours 
of the integrated intensities of the NH 3 (1,1) main hyperfine line. The contours are in steps of 5 mjy beam -1 km s —1 x [3, 6 , 9, 12, 15, 
20, 25]. The dash circles represent the 2' primary beams of the VLA observations. ‘Core 2’ in MM1 is marked with a star. The four 
H 2 O masers are marked with crosses and labelled as W1-W4. (b) Temperature map using the VLA data only, overlaid with contours of 
the integrated intensities of the NH 3 (1,1) main hyperfine line. The contours are in steps of 5 mjybeam -1 kms - 1 x[3, 6 , 9, 12, 15, 20]. 
The three positions exhibiting NH 3 (3,3) emission are marked with circles and labelled as T1-T3, with their spectra and the best gaussian 
fitting shown in the inset. Note that the two maps are plotted using the same color scale from 11 K to 21 K to show their similarity, but 
the highest temperature in (a) is 19 K. 
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Fig. 3.— Two-components fitting of mean spectra of the cores in MM1. The results are listed in Table 3. The NH 3 (1,1) spectra are 
between 70 and 100 kms -1 , while the NH 3 (2,2) spectra are shifted by 30 kms -1 to be between ~105 and 130 kms -1 , so that they can 
be fitted simultaneously. For each spectrum, the dashed magenta and cyan curves represent the two components that are fitted, while the 
solid red curve represent the sum of them. The horizontal solid lines mark the 3cr levels of each spectrum. 



Fig. 4.— Two-components fitting of mean spectra of MMl-pl, using the VLA data only. The legends are the same as in Figure 3. 
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Fig. 5.— The SMA 1.3 mm continuum emission of MM1, MM5 and MM7/8. The dashed circles in each panel mark the 56 ;/ primary beams 
of the SMA. The ellipses mark FWHM of the 2D gaussians fitted to each core, (a) The contours are in steps of 1.1 mjy beam -1 x[—5, —3, 3, 
5, 7, 9, 15, 25, 35, 45]. The two H 2 O masers, W1 and W4, are marked with crosses, (b) The contours are in steps of 0.9 mJy beam -1 x[—3, 
3, 5, 7]. (c) The contours are in steps of 0.81 mjy beam -1 x [—3, 3, 5, 7]. 
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5.5 216.0 216.5 217.0 219.5 220.0 220.5 221.0 225.5 226.0 226.5 227.0 229.5 230.0 230.5 231.0 


O 

MM5-p1 O 
x 



MM7-p1 




MM8-p1 


O 

O 


O 

O 




117.0 217.5 218.0 218.5 219.0 219.5 220.0 220.5 229.0 229.5 230.0 230.5 231.0 231.5 232.0 232.5 


Sky Frequency (GHz) 


lines of MMl-pl, MM5-pl, MM7-pl, and MM8-pl detected by the SMA. The detected species are labelled in each 
1 dashed lines mark the 3cr levels of each spectrum. 
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Fig. 7.— SMA CO channel maps of MM1. The contours are in steps of 0.1 Jy beam -1 x [-5,5,10,15,20,25,30]. MMl-pl and MMl-p2 are 
marked with two stars. The two CO outflows are represented by arrows, with colors indicating blue or red shifted components. The three 
H 2 O masers, Wl, W3, and W4 are marked with asterisks. 



40 20 0 -20 -40 

RA offset (arcsec; J2000) 


Fig. 8.— SMA CO channel maps of MM5. The contours are in steps of 0.05 Jy beam 1 x [-5,5,10,15]. The three cores in MM5 are marked 
with stars. The arrows represent the CO outflow with colors indicating blue or red shifted component. 
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Fig. 9. — Left: Virial parameters of cores in G28.53 as well as in the hot core sample of Lu et al. (2014). Constant radial density profiles 
are assumed for all cores. Right: FWHM linewidths vs. rotation temperatures of cores in G28.53 as well as in the hot core sample of Lu 
et al. (2014). 



18 h 44 m 18.6 s 18.2* 17.8* 18 h 44 m 17.B* 17.4* 17.0* 16.6" 


RA (J2000) RA (J2000) 

Fig. 10.— (a) SMA CH 3 OH velocity field of MMl-pl, overlaid with its integrated intensity. The contours are in steps of 
0.13 Jybeam -1 kms —1 x[3, 5, 7, 9, 11]. The star marks MMl-pl. The red and blue arrows represent the CO outflow associated with 
MMl-pl as in Figure 7. (b) SMA H 2 CO velocity field of MM5, overlaid with its integrated intensity. The contours are in steps of 

0.1 Jybeam -1 kms - 1 x[3, 5, 7, 9, 11, 13, 15, 17]. The three stars mark MM5-pl-p3. The red and blue arrows represent the CO outflow 
in MM5 as in Figure 8 . 
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Fig. 11.— GBT NH 3 (2,2) velocity field of MM1 (the 87 kms 1 component) and MM7/8. The 30 // beam of GBT at K band is shown in 
the lower left corner in each panel. Left: MMl-pl/p2 are marked with stars. Right : MM7-pl and MM8-pl/p2 are marked with stars. 




















